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THE NGC 404 NUCLEUS: STAR CLUSTER AND POSSIBLE INTERMEDIATE MASS BLACK HOLE 
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ABSTRACT 

We examine the nuclear morphology, kinematics, and stellar populations in nearby SO galaxy 
NGC 404 using a combination of adaptive optics assisted near-IR integral-field spectroscopy, opti- 
cal spectroscopy, and HST imaging. These observations enable study of the NGC 404 nucleus at a 
level of detail possible only in the nearest galaxies. The surface brightness profile suggests the presence 
of three components, a bulge, a nuclear star cluster, and a central light excess within the cluster at 
radii <3 pc. These components have distinct kinematics with modest rotation seen in the nuclear star 
cluster and counter-rotation seen in the central excess. Molecular hydrogen emission traces a disk with 
rotation nearly orthogonal to that of the stars. The stellar populations of the three components are 
also distinct, with half of the mass of the nuclear star cluster having ages of ~1 Gyr (perhaps resulting 
from a galaxy merger), while the bulge is dominated by much older stars. Dynamical modeling of the 
stellar kinematics gives a total nuclear star cluster mass of 1.1 x 10 7 M©. Dynamical detection of a 
possible intermediate mass black hole is hindered by uncertainties in the central stellar mass profile. 
Assuming a constant mass-to-light ratio, the stellar dynamical modeling suggests a black hole mass of 
< 1 x 10 5 M©, while the molecular hydrogen gas kinematics are best fit by a black hole with mass of 
4.5^2 x 10 5 Mq. Unresolved and possibly variable dust emission in the near-infrared and AGN-like 
molecular hydrogen emission line ratios do suggest the presence of an accreting black hole in this 
nearby LINER galaxy. 

Subject headings: galaxies:nuclei - galaxies: elliptical and lenticular, cD - galaxies: kinematics and dy- 
namics - galaxies: formation - galaxies: structure - galaxies: individual (NGC 404) 



1. INTRODUCTION 

The centers of galaxies contain both massive black 
holes (MBHs) and nuclear star clusters (NSCs). The 
presence of MBHs has been dynamically measured in 
about 50 galaxies and it appears that most massive 
galaxies probably have a MBH (e.g. Richstone et al. 
1998; Graham 2008). Nuclear star clusters are compact 
{r e ff ~ 5 pc), massive (~ 10 r M ) star clusters found 
at the center of a majority of spirals and lower mass 
ellipticals (Boker et al. 2002; Carollo et al. 2002; Cote 
et al. 2006). Unlike normal star clusters, they have mul- 
tiple stellar populations with a wide range of ages (Long 
et al. 2002; Walcher et al. 2006; Rossa et al. 2006; Seth 
et al. 2006; Siegel et al. 2007). Nuclear star clusters coex- 
ist with MBHs in some galaxies (Filippenko & Ho 2003; 
Seth et al. 2008a; Shields et al. 2008; Graham & Spitler 
2009). However, the nearby galaxies M33 and NGC 205 
have NSCs but no apparent central black hole (Gebhardt 
et al. 2001; Valluri et al. 2005), while some high mass core 
elliptical galaxies have MBHs but lack NSCs (Cote et al. 
2006). 

Occupation fractions and masses for MBHs in lower 
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mass galaxies retain the imprint of the seed black holes 
(BHs) from which they form, information which has been 
lost due to subsequent accretion in higher mass galaxies 
(e.g. Volonteri et al. 2008). However, the presence and 
mass of MBHs in lower mass galaxies is very poorly con- 
strained. Most galaxies are too far away to measure the 
dynamical effect of a BH with mass <10 6 M Q (often re- 
ferred to as intermediate mass black holes; IMBHs) with 
current instrumentation. Thus the presence of IMBHs in 
galaxy centers has only been inferred when AGN activ- 
ity is observed (e.g. Filippenko & Sargent 1989; Greene 
& Ho 2004; Satyapal et al. 2007). These AGN provide 
only a lower limit on the number of MBHs in lower mass 
galaxies, and the BH mass estimates from the AGN are 
quite uncertain. 

Dynamical measurements of MBHs in nearby massive 
galaxies have revealed that the mass of a galaxy's cen- 
tral MBH is correlated with its bulge mass (Kormendy 
& Richstone 1995; Magorrian et al. 1998; Haring & Rix 
2004) . The scaling of BH mass with the large-scale prop- 
erties of galaxies extends to many measurable quantities 
including the bulge velocity dispersion (Ferrarese et al. 
2000; Gebhardt et al. 2000; Graham 2008; Giiltekin et al. 
2009). More recently, Ferrarese et al. (2006), Wehner & 
Harris (2006), and Rossa et al. (2006) have presented ev- 
idence that NSCs scale with bulge mass and dispersion in 
elliptical and early-type spiral galaxies in almost exactly 
the same way as MBHs. 

The similarity between the NSC and MBH scaling rela- 
tionships led Ferrarese et al. (2006) and Wehner & Harris 
(2006) to suggest that MBHs and NSCs are two differ- 
ent types of central massive objects (CMO) both of which 
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contain a small fraction (~0.2%) of the total galaxy mass. 
This correlation of the CMO mass with the large scale 
properties of galaxies suggests a link between the forma- 
tion of the two. However, the nature of this connection 
is unknown, as is the relation between NSCs and MBHs. 
Studies of NSCs can help address these issues. Their 
morphology, kinematics, and stellar populations contain 
important clues about their formation and the accretion 
of material into the center of galaxies (e.g. Hopkins & 
Quataert 2010). 

The scaling relationships of large-scale galaxy prop- 
erties with the mass of the CMO might indicate that 
galaxy centers should be simple systems. The Milky 
Way center clearly shows this not to be the case. It is 
an incredibly complicated environment with a black hole 
(M BH = 4 x 10 6 M Q ; Chez et al. 2008) and a nuclear 
star cluster {M NSC ~ 3 x 10 7 M© Genzel et al. 1996; 
Schodel et al. 2007; Trippe et al. 2008) that contains 
stars of many ages and in different substructures includ- 
ing disks of young stars in the immediate vicinity of the 
black hole (e.g. Lu et al. 2009; Bartko et al. 2009). Only 
through understanding these complex structures in the 
Milky Way and finding other examples in nearby galaxies 
can we hope to fully understand MBH and NSC forma- 
tion, and the links between these objects and their host 
galaxies. 

In our current survey, we are looking at a sample of the 
nearest galaxies (D < 5 Mpc) that host NSCs. We are 
resolving their properties using a wide range of observa- 
tional data to understand how NSCs form and their rela- 
tion to MBHs. In our first paper on nearby edge-on spiral 
NGC 4244 (Seth et al. 2008b), we showed evidence that 
the NSC kinematics are dominated by rotation, suggest- 
ing that it was formed by episodic accretion of material 
from the galaxy disk. 

This paper focuses on the NSC and possible MBH in 
NGC 404, the nearest SO galaxy. Table 1 summarizes 
its properties. The stellar populations of the galaxy are 
predominantly old and can be traced out to 600" (9 kpc) 
(Tikhonov et al. 2003; Williams et al. 2010). However, 
HI observations show a prominent nearly face-on HI disk 
at radii between 100-400", with detectable HI out to 
800" (del Rio et al. 2004). CO observations and optical 
color maps show the presence of molecular gas and dust 
within the central ~20" (300 pc) of the galaxy (Wiklind 
& Henkel 1990; Tikhonov et al. 2003), lying primarily to 
the NE of the nucleus. A nuclear star cluster in the cen- 
tral arcsecond of NGC 404 was noted by Ravindranath 
et al. (2001) from an analysis of NICMOS data. 

The presence of an AGN in NGC 404 is a controver- 
sial topic. The optical spectrum of the nucleus has line 
ratios with a LINER classification (Ho et al. 1997). As- 
suming a distance of ^3 Mpc, NGC 404 is the nearest 
LINER galaxy; other nearby examples include M81 and 
NGC 4736 (M94). Recent studies have shown that most 
LINERS do in fact appear to be AGN, with a major- 
ity of them having detected X-ray cores (Dudik et al. 
2005; Gonzalez-Martin et al. 2006; Zhang et al. 2009), 
radio cores (Nagar et al. 2005), and many of them poss- 
csing mid-IR coronal lines (Satyapal et al. 2004) and UV 
variable cores (Maoz et al. 2005). However, NGC 404 
is quite unusual in its properties. No radio core is ob- 
served down to a limiting flux of 1.3 mJy at 15 GHz 



TABLE 1 
NGC 404 Properties 



Distance" 


3.06 Mpc 


TO — M 


27.43 


pc/" 


14.8 


Galaxy M vfi , M I}H b 


-17.35, -18.36 


Bulge/Total Luminosity* 2 


0.76 


Bulge M V}0 , M Ifi d 


-17.05, -18.06 


Bulge Mass 6 


9.2x10 s M 


HI Gas Mass-f 


1.5x10 s M 


Molecular Gas Mass 9 


6xl0 6 M 


Central Velocity' 1 


-58.9 kms" 1 



Note. — (a) TRGB measurement from Karachentsev et al. 
(2004), (b) at r < 200" from Tikhonov et al. (2003) corrected 
for foreground reddening, (e) for r < 200" from Baggctt et al. 
(1998), (d) combining above values; we derive a slightly brighter 
bulge M I}0 = -18.20 in §3, (e) using M/Lj = 1.28 from §5.4, (/) 
del Rio et al. (2004), (g) Wiklind & Henkel (1990), corrected to 
D = 3.06 Mpc, (h) heliocentric velocity, see §4.1. 

with the VLA in A array (Nagar et al. 2005), however, 
del Rio et al. (2004) do detect an unresolved 3 mJy con- 
tinuum source at 1.4 GHz using the C array. A com- 
pact X-ray source is detected (Lira et al. 2000; Eracleous 
et al. 2002), but its low luminosity and soft thermal spec- 
trum indicates that it could be the result of a starburst. 
Signatures of O stars are seen in the UV spectrum of 
the nucleus, however dilution of these lines suggests that 
^60% of the UV flux could result from a non-thermal 
source (Maoz et al. 1998). This suggestion is supported 
by more recent UV observations which show that the 
UV emission is variable, declining by a factor of ~3 be- 
tween 1993 and 2002 (Maoz et al. 2005). HST observa- 
tions of Ha show that the emission occurs primarily in a 
compact source 0.16" north of the nucleus and in wispy 
structures suggestive of supernova remnants (Pogge et al. 
2000). The [O III] emission has a double lobed struc- 
ture along the galaxy major-axis (Plana et al. 1998), and 
has a higher velocity dispersion than Ha near the galaxy 
center (Bouchard et al. 2010). The mid-IR spectrum of 
NGC 404 shows evidence for high excitation consistent 
with other AGN (Satyapal et al. 2004). In particular 
the ratio of the [O IV] flux relative to other emission 
lines ([Ne II], [Si II]) is higher than any other LINERs in 
the Satyapal et al. (2004) sample and is similar to other 
known AGN. However, [NeV] lines, which are a more re- 
liable indicator of AGN activity (Abel & Satyapal 2008), 
are not detected. In summary, the case for an accreting 
MBH in NGC 404 remains ambiguous, with the variable 
UV emission providing the strongest evidence in favor of 
its existence. 

In this paper we take a detailed look at the central 
arcseconds of NGC 404 and find evidence that it contains 
both a massive nuclear star cluster and a black hole. 
In §2 we describe the data used in this paper including 
adaptive optics Gemini/NIFS observations. We then use 
this data to determine the morphology in §3, the stellar 
and gas kinematics in §4, and the stellar populations in 
§5. In §6 we present dynamical modeling from which we 
derive a NSC mass of 1.1 x 10 7 M Q and find mixed results 
on detecting a possible MBH with mass < 10 6 M©. In 
§7 we discuss these results, concluding in §8. 
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2. OBSERVATIONS & REDUCTION 

This paper presents a combination of a wide variety of 
data including: (1) high spatial resolution near- infrared 
IFU spectroscopy with Gemini NIFS, (2) optical long- 
slit spectra from the MMT, and (3) HST imaging data 
from the UV through the NIR. Using these data we have 
determined the morphology, gas and stellar kinematics, 
and stellar populations of the central regions of NGC404. 
In this section we describe these data and their reduction. 

2.1. Gemini NIFS data 

Observations of NGC 404 with Gemini's Near-Infrared 
Integral Field Spectrometer (NIFS) are part of an on- 
going survey of nearby nuclear star clusters. The data 
were taken using the Altair laser-guide star (LGS) sys- 
tem on the Gemini North telescope on the nights of 21st 
and 22nd September, 2008. The nucleus itself was used 
as a tip-tilt star. A total of twelve 760 second exposures 
were obtained; half were dithered on-source exposures 
while the others were sky exposures offset by 131" to 
the northwest taken after each on-source exposure. A 
telluric calibrator, HIP 1123, an Al star, was observed 
both nights at an airmass similar to the observations. 

The data were reduced using the Gemini v. 1.9 IRAF 
package, utilizing pipelines based on the NIFSEXAM- 
PLES scripts. Ronchi mask and arc lamp images were 
used to determine the spatial and spectral geometry of 
the images. The sky images were subtracted from their 
neighboring on-source images, and then these images 
were flat-fielded, had bad pixels removed, and were split 
apart into long-slit slices. The spectra were then di- 
vided by the telluric calibrator (using NFTELLURIC) 
after correcting for the Br7 absorption in the calibra- 
tor. We made minor alterations to the Gemini pipeline 
to enable propogation of the variance spectrum. This 
also required creation of our own IDL version of the 
NIFCUBE to make the final data cubes. In this pro- 
cess the initial pixel size of the data, 0.103" in the verti- 
cal direction and 0.045" in the horizontal direction, was 
rcbinned to 0.05" x 0.05" pixels in the final data cube 
for each on-source exposure. Data cubes for the six 
dithered exposures were then shifted spatially and com- 
bined with cosmic-ray rejection using our own IDL pro- 
gram. Data from the 2nd night was spectrally shifted 
to the wavelength solution of the first night before com- 
bining. The heliocentric correction was nearly identical 
on both nights (+13.25 kms" 1 ). The NIFS spectrum 
within a 1" aperture is shown in Fig. 1 and shows strong 
absorption and emission lines. 

We also processed the sky frames in a manner iden- 
tical to the on-source exposures. We used the resulting 
data to calibrate the spectral resolution of our observa- 
tions. Using thirteen isolated, strong, doublet sky lines, 
we measured the spectral resolution of each pixel in our 
final data cube. We found significant variations in the 
spectral resolution across the detector, with a median 
resolution of 4.36 A FWHM (A/ A A = 5275), and values 
ranging from 3.9 to 5.1 A FWHM. 

A if-band image was made by collapsing all chan- 
nels after multiplying by (1) a 9500K black-body to 
correct for the shape of the telluric spectrum, and (2) 
the 2MASS response curve. Flux calibration was ob- 
tained from 2MASS by deriving a zeropoint from im- 



ages of both sources and calibrators in our sample. We 
flux calibrated our spectra using synthetic 2MASS pho- 
tometry of our telluric calibrator HIP 1123, assuming 
its spectral shape resembles that of a 9500K black body. 
From this flux calibration we obtained a flux calibrated 
spectrum of our NGC404 NIFS data cube with units of 
ergs s _1 cm~ 2 A -1 . The flux calibration of the image and 
spectrum is good to ~10%. 

2.1.1. PSF Determination 

Understanding the point spread function (PSF) is im- 
portant for interpreting the kinematic observations pre- 
sented in §4, and is even more critical for the dynamical 
modeling in §6. We used two methods for measuring the 
PSF of the calibrator: (1) we created images from the 
telluric calibrator for both nights to estimate the PSF, 
and (2) we compared the final NIFS if-band images of 
the NGC404 nucleus to HST F814W data. 

We initially tried modeling the PSF as a double Gaus- 
sian (as in Krajnovic et al. 2009). However, double Gaus- 
sians were a very poor fit to the wings of the PSF in the 
telluric calibrator images. Using an inner Gaussian + 
outer Moffat profile (£(r) = E /[(l + (r/r d ) 2 } 4 - 765 ) ad- 
equately described the telluric calibrator PSF. The best 
fit on both nights was a PSF with an inner Gaussian 
FWHM of 0"12, and an outer Moffat profile r d of 0'.'95, 
each containing about half the light. The two nights 
PSFs were very similar (< 10% differences) in both shape 
and fraction of light in each component. Unfortunately, 
our telluric images were not well-dithered and thus the 
central PSF is not well sampled. 

We further refined the PSF by fitting a convolved ver- 
sion of the F814W band HST image to the NIFS K- 
band image. Although HST/NIC2 and NIC3 F160W 
data were also available, the larger pixel scale and the 
limited resolution of HST in the NIR (FWHM-0.14") 
make the F814W data a better image to use as a PSF 
reference. Due to the dust around the nucleus, we fit 
only the western, apparently dust-free, half of the nu- 
cleus (Fig. 2). Also, we find in §3.1 that the NIFS image 
contains compact non-stellar emission from hot dust at 
the center; we therefore correct the if-band image for 
this dust correction before determination of the PSF. Fi- 
nally, because the profile of the galaxy is much shallower 
than the PSF, we fix the functional form of the outer pro- 
file to that derived for the telluric star, but allow the in- 
ner PSF width and fraction of light in each component to 
vary. The best-fitting PSF has a core with FWHM- 0'/09 
(assuming an HST resolution of 0"065 based on TinyTim 
models; Krist (1995)) containing 51% of the light, with 
the outer Moffat (with r d = 0"95) function containing 
49% of the light. 

This PSF suffers from a number of uncertainties. Most 
notably, the presence of a bluer population at the center 
of the cluster implied by Fig. 2 may make the F814W 
image more centrally concentrated relative to the dust- 
corrected if-band image. Due to our fitting the NIFS 
PSF core to a convolution of the F814W image, this effect 
could lead to an over-estimate in the core FWHM. 

2.1.2. Stellar Kinematics Determination 

Derivation of the stellar kinematics was done as de- 
scribed in Scth et al. (2008b) with a couple of additions. 
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Fig. 1. — Left - Optical spectrum of the NGC 404 nuclear star cluster (l" X 2. 3" aperture) taken with the blue channel spectrograph on 
the MMT. Right - fC-band spectrum of the NGC 404 nuclear star cluster using a 1" radius aperture from our Gemini/NIFS data. Emission 
lines of H2 are labeled above the spectrum (see §3.2), while strong absorption features are labeled below the spectrum (including Br7 which 
is seen weakly in absorption and emission, see Fig. 6). 



In brief, we first spatially binned the data to a target 
signal-to- noise (S/N) of 25 using the Voronoi tessela- 
tion method described in Cappellari & Copin (2003). 
The S/N of the data was as high as 175 in the cen- 
ter, and data within ~ Of! 7 was left unbinned. Next, 
we determined the stellar line-of-sight velocity distribu- 
tion (LOSVD) using the penalized pixel-fitting of Cap- 
pellari & Emsellem (2004) to derive the velocity, disper- 
sion, h3, and M components in the wavelength region 
from 22850A to 23900A. We used high-resolution tem- 
plates from Wallace & Hinkle (1996) which has super- 
giant, giant and main sequence stars between spectral 
types G and M; we used eight of these stars with the 
most complete wavelength coverage. These templates 
were convolved to the instrumental resolution measured 
in each pixel from the sky lines before fitting the LOSVD. 
Nearly identical kinematic results and fit qualities were 
obtained using the larger GNIRS library (Winge et al. 
2008) as templates. We performed LOSVD fits over dif- 
ferent ranges of wavelengths and with different subsets 
of templates and achieved consistent results, suggesting 
template mismatch is not a significant issue. The higher 
spectral resolution and more accurate velocity zeropoints 
made the Wallace & Hinkle (1996) better than the Winge 
et al. (2008) templates for deriving the kinematics. Er- 
rors on the LOSVD were calculated by taking the stan- 
dard deviation of values derived from 100 Monte Carlo 
simulations in which appropriate levels of Gaussian ran- 
dom noise was added to each spectral pixel before remea- 
suring the LOSVD. Errors on the radial velocities ranged 
from 0.5 to 6 kms" 1 depending primarily on the S/N. 

2.2. Optical Spectroscopy 

Optical spectroscopy of the NGC 404 cluster was ob- 
tained (along with spectra of other nearby NSCs) with 
the blue channel spectrograph (BCS) on the MMT 6.5m 
telescope on the night of 19 Jan 2009. We used the 500 



line grating and a 1 arcsecond slit to obtain a resolution 
of ~3.6A FWHM between 3680A and 6815A with a pixel 
size of 0'.'3 spatially and 1.17A spectrally. Two 1200s ex- 
posures were taken, with one being affected by clouds. 
Due to variations in the slit width, flat-fielding of the 
data was accomplished by carefully constructing a com- 
posite flat field using the pattern from quartz flats taken 
just after the exposure with pixel-to-pixel information 
drawn from flats throughout the night. Wavelength cal- 
ibration was done using arcs taken just after the science 
exposures. Spectral extraction was accomplished using 
the standard IRAF DDSLIT routine using optimal extrac- 
tion over the central 2'.' 3 of the slit while sky regions were 
selected to avoid regions of emission near the nucleus. 
Guide camera seeing was around 0'.'8, however within the 
instrument it appeared to be significantly larger (~ l'/7). 
Three flux calibrators, Gliese 191B2B, HZ44, and VMa2 
(Bessell 1999) were observed on the same night with the 
same observing set up. Flux calibration was achieved 
by combining information from all three flux calibrators 
using SENSFUNC. Final residuals between the fitted sensi- 
tivity function and the 3 calibrators were <0.03 mags be- 
tween 3680A and 6200A. We note that 2nd order contam- 
ination was present beyond 6200A. One of the NGC 404 
exposures and 2/3 calibrator exposures were affected by 
clouds. However, by scaling to the clear exposures, the 
absolute calibration of the spectra appears to be quite 
good; comparison of emission line fluxes to those in Ho 
et al. (1997) shows agreement to within 5% through a 
similar size aperture. The S/N per pixel of the final 
combined spectrum rang ed from 100 at 3700A to 250 
around 5000A. The extracted nuclear spectrum is shown 
in Fig. 1. 

2.3. HST data 

A wide variety of HST data are available for the nu- 
cleus of NGC 404 taken with the FOS, WFPC2, NIC- 
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Fig. 2.— The HST WFPC2 F547M - F8UW (approximately 
V — I) color map of the nuclear region. Contours show the 
NIFS continuum emission at fi K of 11.2, 12.2, 13.2, 14.2 and 15.2 
mag/arcsec 2 . The image is centered on the dust-corrected peak of 
the continuum light with North up and East to the left. Redder 
regions result from dust extinction. The western half of the nucleus 
appears to have little internal extinction (see §5.2). 

MOS and ACS HRC and WFC. In this paper we will 
use WFPC2 data in the F547M, F656N (PID: 6871), 
F555W, and F814W (PID: 5999) filters, ACS WFC data 
in F814W filter (PID: 9293) and NICMOS NIC2 and 
NIC3 data in the F160W filter (PID: 7330 & 7919). Data 
were downloaded from the HST archive and the Hubble 
Legacy Archive (HLA). The WFPC2 data all had the 
nucleus in the PC chip and were unsaturated, while the 
ACS F814W data is saturated at the center but covers 
the galaxy to larger radii than the WFPC2 data. Astro- 
metric alignment of this data and the NIFS data was tied 
to the F814W WFPC2 data. The centroided position of 
the nucleus was used to align each image to the F814W 
data. Although dust clearly affects the area around the 
nucleus, the F547M-F814W color map (Fig. 2) and the 
alignment of sources in the NIFS Br7 map with the HST 
Ha image (see §3.2) suggests that the very center of the 
galaxy does not suffer from significant internal dust ex- 
tinction. 

3. MORPHOLOGY 

To determine the extent and luminosity of the nuclear 
star cluster in NGC 404, it is necessary to fit the rela- 
tively steep surface brightness profile of the inner part of 
the galaxy. Using Sersic fits to I and H/K-b&nd data, 
we find evidence that a nuclear star cluster dominates 
the central 1" (15 pc) of the galaxy, and that an addi- 
tional emission component appears to be present within 
the central few parsecs (see Fig. 3). 

Before describing our results on the central parts of the 
galaxy, we summarize previous studies of the large scale 
light-profile of NGC 404. Using ground-based V-band 
data out to nearly 300", Baggett et al. (1998) decompose 
the SB profile of NGC404 into a r 1 / 4 bulge and exponen- 
tial disk. They find a bulge effective radius of 64" and a 
disk with scale length of 130". The bulge dominates the 
light out to ^150"and the bulge-to-total ratio is ^0.7. 
Tikhonov et al. (2003) found that the number counts of 
red giant branch stars between 100-500" (1.5-7.5 kpc) are 
also well explained by an exponential disk with a similar, 



^106" scale length. 

We used archival HST data, our NIFS if -band image, 
and the 2MASS Large Galaxy Atlas (Jarrett et al. 2003) 
to construct the two surface brightness profiles shown in 
Fig. 3. Surface brightness measurements were obtained 
on all images using the ellipse program in IRAF. The 
F814W-band (I) profile combines HST WFPC2 data 
from the PC chip at radii below 5" with ACS/WFC 
data at larger radii (the ACS image was saturated at 
smaller radii). Diffraction spikes from /3And (7' to the 
SE) and bright foreground stars and background galax- 
ies were masked. Dust is obviously present near the cen- 
ter and this was masked out as well as possible using a 
FhhhW - F814W color map at radii between 0"3 and 
5" from the center. Comparing this profile with the NIR 
light profile suggests that some dust extinction remains 
within the inner 1". The NIR profile combines the NIFS 
if-band image (at r < 1"), with the F160W-band (H) 
NICMOS/NIC3 image (1 < r < 5") and H band data 
from the 2MASS Large Galaxy Atlas (at r > 5"). The 
NIFS data was scaled to match the fZ-band zeropoint; no 
scaling was needed to match the 2MASS and NICMOS 
data. The final F814W and NIR surface brightness pro- 
files are shown in Fig. 3. 

We first fit the outer part of both profiles to a Sersic 
profile (see Graham & Driver 2005, for a complete discus- 
sion of Sersic profiles). During our fitting, the Sersic pro- 
files were convolved with the PSF. The 2MASS data was 
only usable to 70" due to low S/N, and even within this 
radius the sky subtraction is suspect due to nearby /3And. 
The F814W data is more reliable at these larger radii, 
however we fit only the measurements within 80" due to 
uncertainty in the sky background and possible contribu- 
tion from a disk at larger radii. At radii beyond 5", the 
profiles are well fit by similar Sersic profiles with n ~ 2.5 
and effective (half-light) radius of ^43" (640 pc). Exact 
values for all fits are given in Table 2. We identify this 
outer component as the bulge given the results of Baggett 
et al. (1998), who show the bulge component dominates 
at r < 150". It was not possible to fit either profile down 
to smaller radii with a single Sersic. 

Following Graham & Guzman (2003); Cote et al. 
(2006), we identify the excess light over the bulge at small 
radii to be a NSC. We discuss the nature of this compo- 
nent more in §7.1. In both profiles, the NSC component 
is brighter than the bulge at r < 1". We fit this com- 
ponent with an additional Sersic component (following 
Graham & Spitler 2009) and find it has an effective ra- 
dius of - 0"7 (10 pc). The NSC Sersic n value is quite 
sensitive to the functional form of the PSF, therefore we 
prefer the F814W value of n=1.99 due to the better- 
known PSF. These parameters are in the ranges typical 
of NSCs (Cote et al. 2006; Graham & Spitler 2009) in 
other galaxies. The absolute magnitude of this compo- 
nent (based on the model) is Mr, Mh = -13.64, -15.34 
within the central 3", and M r ,M H = -13.17,-14.93 
within 1". 

In the inner 0"2, there is a clear excess luminosity 
above the double-Sersic fit in both profiles. Adding up 
the light in excess of the profile we find an Mj , Mh = 
— 11.81,-12.61. Note that these magnitudes are quite 
sensitive to the PSF and Sersic n indices of the NSC 
fit. We will show below that in the NIR, a fraction of 
this light appears to be emission by hot 950 K dust. A 
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Fig. 3. — Surface brightness profiles of NGC 404 created from HST, Gemini/NIFS, and 2MASS data. In both panels, the observed surface 
brightness profiles are plotted with symbols, the solid line shows the best-fit two-Sersic model, the dashed line the outer (bulge) component, 
and the dot-dashed line the inner (NSC) component. The PSF is shown as a dotted red line in both plots. Residuals disappear off the 
plot at small radii to the unmodeled extra-light component. Left - the F8f4W profile constructed from both WFPC2/PC (diamonds) and 
ACS/WFC (squares) data. This data is the most reliable data available at large radii. Right - the NIR profile constructed from three 
images: NIFS X-band (diamonds), HST NICMOS/NIC3 (triangles) and 2MASS large galaxy atlas data (squares Jarrett et al. 2003). Blue 
diamonds show a correction for hot dust emission discussed in §3.1. 



correction for this dust emission is shown with blue di- 
amonds in Fig. 3. However, dust alone cannot explain 
the luminosity of this excess light, accounting for only 
^35% of the excess in the NIR profile. Furthermore, at 
/-band, no significant dust emission would be expected 
at all. Two additional clues paint a complicated picture. 
First, the color in the central 0"2 is somewhat bluer in 
the HST F547M-F8UW images (Fig. 2) and has com- 
pact UV emission (Maoz et al. 2005), both suggestive of 
the presence of young stars. Second, the stars in this 
area appear to be counter-rotating relative to the NSC 
and galaxy (see §4.1). Thus the central excess appears to 
result from a combination of hot dust emission, younger 
stars, and perhaps AGN continuum emission in the op- 
tical and UV. 

In summary, we find evidence for 3 components at the 
center of NGC 404. 

1. A bulge with Sersic n ~ 2.5 and effective radius 
of 640 pc that dominates the light between 1" and 
80" (and beyond). 

2. A nuclear star cluster which dominates the light in 
the central arcsecond (15 pc). 

3. A central excess at r < 0'.'2. This excess appears 
to result from the combination of hot dust, young 
stars, and perhaps AGN continuum emission. 

Within the central 1" in the /-band (roughly the area 
included in the MMT spectrum analyzed in §5), ^55% of 
the light comes from the NSC, 29% from the underlying 
bulge, and 16% from the central excess. 

The isophotes in the ACS F814W data at 10-80"have a 
position angle of ~80 degrees and an ellipticity of ~0.08. 
This matches well with the morphology of the HI disk on 
a large scale; del Rio et al. (2004) finds an HI inclination 
of ~10 degrees and a position angle of about ~80 degrees 



TABLE 2 
Morphological Fits 



F814 W Profile 

Component Sersic n Teft [''I Me f f,FSi4,w [ ma g/" 2 ] Afj 

Bulge 2A3 CU> 20718 -18.09 

NSC 1.99 0.74 16.02 -13.64 

Central Excess -11.81 

NIR Profile 

Component Sersic n r e tf ["] l^effH [ ma g/" 2 ] Mh 

■BUTgi 2T7i itr ' 18.82 ^mwi 

NSC 2.61 0.68 14.15 -15.34 

Central Excess -12.61 
Central Excess (Dust Corrected) -12.14 



Note. — Bulge magnitudes given for r < 80", and NSC magni- 
tudes for r < 3". Central excess magnitudes are for the total flux 
above the fit in the central arcsecond. Distance modulus is assumed 
to be 27.43 giving 1"=14.8 pc. Note that the magnitudes are not 
corrected for foreground extinction (At = 0.11, Ah = 0.03; Schlegel 
et al. 1998). 

for radii between 100-300". Due to dust absorption, the 
clearest picture of the inner arcsecond of NGC 404 comes 
from the NIFS ZT-band image. As can be seen from the 
contours in Fig. 2, the position angle at r < 0"25 has a 
PA of 60-70 degrees. The PA then abruptly changes to 
25 degrees at r — 0"3 and gradually increases to the large 
scale PA of 80 degrees by ~1". The ellipticity decreases 
from ~0.2 near the center to 0.07 at 1". These changes 
fit well with our picture of three distinct components, 
each having their own distinct morphology. 

3.1. A Hot Dust Component 

Within the central 0.2" , the NIFS spectra show a sig- 
nificant spectral flattening and a reduction in the depths 
of the CO lines by about 10%. A comparison of the cen- 
tral spectrum to single stellar populations models from 
Maraston (2005) shows that its spectral slope is flat- 
ter/redder than any stellar population models. Further- 
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Fig. 4. — Fit of the central NIFS spectrum of NGC 404 to a model 
including emission from hot dust. The black line shows the central 
spectrum, with the overplotted line showing the best fit. This fit 
consists of a template stellar spectrum taken from an annulus of 
radius O'.'5-l'.'O, plus a black-body of temperature 950 K (shown 
in gray). Note that emission lines were excluded from the fit and 
plot. 

more, the younger stellar populations expected in the 
central excess should have a somewhat bluer, not redder, 
continuum. This redder continuum suggests two possi- 
bilities (1) strong dust absorption or (2) emission from 
hot dust. 

We rule out dust absorption as a possibility because 
the required dust absorption to flatten the spectrum is 
Ay ~ 6 magnitudes assuming the reddening law from 
Cardelli et al. (1989). The F547M- F8UW shows that 
the nucleus is actually bluer than the surrounding areas 
(F547M - F8UW ~ 0.8). Furthermore, the area to the 
west of the nucleus appears unreddened with F547W — 
F814W ~ 0.95, consistent with the expected color of 
the stellar population based on spectroscopic fits (see 
§5.2), while the eastern half of the nucleus is at most 0.8 
magnitudes redder in F5A7M — F81AW (approximately 
V— I). This suggests a maximum extinction in the region 
around the nucleus of Ay ~ 2. 

Therefore, hot dust appears to be the cause of the flat- 
tened continuum in the central pixels. Hot dust is a com- 
mon feature in Seyfert nuclei at NIR wavelengths, and 
is typically found to have temperatures of 800-1300 K 
(Alonso-Herrero et al. 1996; Winge et al. 2000; Riffel 
et al. 2009). Studies of Seyfert galaxies with NICMOS by 
Quillen et al. (2000, 2001) show that the hot-dust com- 
ponent is frequently unresolved and variable on months 
to year timescales. The presence of hot dust components 
in LINER galaxies appears to be unstudied. 

We fit the contribution and temperature of this hot 
dust emission by assuming that the central spectrum is 
made up of a nuclear star cluster spectrum drawn from 
an annulus of O'.'5-l'.'O plus a pure blackbody spectrum. 
As shown in Fig. 4, the central spectrum (r < 0'.'05) 
is best fit by a dust temperature of 950 K, contribut- 
ing about 20% of the light. The residuals between the 
central spectra and model have a standard deviation of 
1.4% across the band and the fit completely eliminates 
the difference in spectral slopes between the annular and 
central spectrum. 

Fitting spectrum in radial annuli, we find no evidence 
for dust emission beyond 0'.'3 from the nucleus. Assum- 
ing a uniform dust temperature of 950 K, we calculated 
the contribution of dust in each individual spaxel. Using 
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Fig. 5.— The surface brightness profile of the NICMOS F160W 
observations (NIC2 in red, NIC3 in orange), 2MASS ii-band 
(dashed), and NIFS if -band observations scaled to match the 
F160W NICMOS observations at r > 1". The dashed line shows 
the NIFS observations after correction for dust emission. 

this, we correct the -ftf-band image for the hot dust con- 
tribution. We use the dust- corrected K-band image and 
central position for all our subsequent analysis. The dust 
emission accounts for only ^35% of the central excess in 
the NIR, as is shown with the blue points in right panel 
of Fig. 3. 

The dust emission is slightly offset to the NE of the 
photocenter. This offset is in the direction of the dust 
absorption clearly seen on larger scales to the E of the 
cluster in the HST color image (Fig. 2). The corrected 
image has a photocenter offset by about a half-pixel from 
the uncorrected image, while the dust emission is off- 
set by ~1.5 pixels from the corrected photocenter. The 
small shift in the photocenter reduces the asymmetry in 
the velocity field, moving the photocenter to the center 
of the central counter-rotation (see §4.1). As found in 
other AGN, the hot dust component is unresolved, with 
a measured FWHM of 0'.'14 (single-Gaussian PSF fits 
give a FWHM of - 0'.'18 for our NIFS data). 

The total dust emission in K band has a magni- 
tude of K ~ 15.4 which translates to a flux of 1.3 x 
10~ 13 crgss -1 cm~ 2 and luminosity of 1.4x 10 38 ergss -1 . 
Assuming the hot dust emission is the result of accretion 
onto a BH of ~ 10 5 M Q (§6), this suggests Lboi/LEdd ^ 
10~ 5 ; typical LINERS have L hol jL Edd - 3 x 10~ 5 (Ho 
2004). The hot dust luminosity is somewhat larger 
than the [O III] luminosity (assuming D=3.06 Mpc) 
of 2 x 10 37 ergss- 1 (Ho et al. 1997) and the possi- 
ble hard X-ray component with 2-10 keV luminosity of 
- 3 x 10 37 ergs sT 1 noted by Terashima et al. (2002). The 
1.6 /im dust emission, [O III], and hard X-ray luminosi- 
ties of Seyfert galaxies are discussed in detail by Quillen 
et al. (2001). Although the luminosities presented here 
are much lower than those in their sample (as would be 
expected for a system with a lower accretion rate), the 
relative luminosities in these three bands fall within the 
scatter of their data suggesting a comparable SED. 

We also find tentative evidence for variability in the 
dust emission between our NIFS images taken in 2008 
and HST NICMOS images taken in 1998. This variabil- 
ity may be expected given the variability in the nuclear 
UV flux reported by Maoz et al. (2005). They find that 
the flux at 2500A decreasing by a factor of ~3 between 
1994 and 2002. However, this comparison is based on 
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Fig. 6. — Le/f - a WFPC2 F656N-F814W map of the nuclear region. Negative values are areas with strong Ha emission. Right - Htj 
emission line map. Negative values mean strong 3vy emission, while positive values result from Br7 absorption. Contours and orientation 
are identical to those in Fig. 2. 



observations by three different cameras, the earliest of 
which was an 0'.'86 aperture spectrum, thus spatial reso- 
lution of the UV variability is not possible. 

Figure 5 shows a comparison of the surface brightness 
profiles of all the NIR data, with the NIFS i^-band image 
scaled to fit the SB profile of the NICMOS F160W im- 
ages at r > 1". Within the central Of! 2, the scaled NIFS 
luminosity is clearly fainter than the NIC2 image. This 
change in brightness could be explained by a bluer stel- 
lar population within the central excess, but the effect is 
much larger than expected - the F547 AI — F814W image 
shows a change of A(F547A/ - F8UW) ~ 0.1 mags be- 
tween the center and the NSC at larger radii, implying a 
A(H-K) < 0.1. However, the observed A(H-K) = 0.4 
after attempting to correct both the NIC2 image and 
NIFS dust-subtracted image for the effects of the PSF. 
Using the original NIFS image (without dust subtrac- 
tion), we still get A(H - K) = 0.3. We suggest that 
this large A(H — K) may result from variable dust emis- 
sion, with the NIC2 observations being taken at an ear- 
lier epoch when the dust-emission was more prominent. 
Multi-epoch HST observations of the nucleus in the UV 
and NIR using the same instrument will be extremely 
useful for verifying a possible variable AGN contribution 
and constraining its SED. 

3.2. NIR Emission Line Excitation 

In this section we discuss the emission lines observed in 
the Gemini/NIFS if-band spectrum. We find that they 
imply that the emission originates in dense molecular gas 
with line ratios typical of those observed in known AGN. 

To obtain estimates for line emission in the NGC 404 
nucleus, we summed together flux-calibrated integrated 
spectra from our NIFS data with radii ranging between 
1 and 25 pixels. We then used both the Wallace & Hin- 
kle (1996) and Winge et al. (2008) templates to fit the 
stellar continuum from 20200A to 24100A excluding the 
parts of the spectrum near known emission lines. Sub- 
tracting this continuum off left us with a pure emis- 



sion spectrum. We measured numerous H2 emission 
lines as well as Br7 emission. The only observable coro- 
nal line commonly seen in AGN is [CaVIII], which falls 
directly on a CO bandhead; an 3 — a upper limit of 
3.6 x 10 -17 ergss -1 cm~ 2 was placed on this line. 

Detected H2 lines (and their vacuum wavelengths) in- 
cluded 1-0 S(2) (20337.6A), 2-1 S(3) (20734.7A), 1-0 S(l) 
(21218.3A), 2-1 S(2) (21542. lA), 1-0 S(0) (22232. 9A), 2- 
1 S(l) (22477.2A) and 1-0 Q(l) (24065.9A). Line ratios 
of these lines can be used to determine the excitation 
mechanisms of the gas. In particular, the ratio of the 
2-1 S(l) and 1-0 S(l) lines is between 0.11 and 0.16 with 
smaller values at larger radii. These values are indicative 
of thermal excitation from shocks or X-ray photons in a 
dense cloud (n > 10 4 cm -3 ), while fluorescent emission 
in a low density medium has a higher 2-1 S(l)/l-0 S(l) of 
^0.55 (e.g. Mouri 1994). The observed ratio is typical of 
those observed in other AGN (e.g. Reunanen et al. 2002; 
Riffcl et al. 2008; Storchi-Bergmann et al. 2009). 

The rotational and vibrational temperatures can be 
estimated from ratios of H2 lines; using the relation 
from Reunanen et al. (2002) T mb = 5600/ln(F(l - 
0S(l))/F(2-lS{l))*1.355) 7 we&ndT mb = 2340±160K 
within an 1" radius aperture. Similarly, the rotational 
temperature can be determined using the relation T rot = 
-1113/ln(F(l-0S(2))/F(l-0S{0))*0.323) (Reunanen 
et al. 2002) which gives a somewhat lower temperature 
of T rot = 1640 ± 360 K within a 1" aperture. The higher 
vibrational temperature may be indicative of some con- 
tribution of fluorescence (e.g. Riffel et al. 2008). 

Unfortunately the Br7 emission line flux can not be ac- 
curately estimated due to uncorrected telluric absorption 
in the vicinity of the Br7 line. However, the distribution 
of Br7 flux is much patchier than the H2 (see Fig. 6) 
suggesting different sites for the two types of emission. 
Note that in starbursting galaxies, the Br7 is stronger 
than the H 2 1-0 S(l) line (Rodriguez- Ardila et al. 2004), 
a possibility which is clearly excluded by our data (see 
Fig. 1). This again is consistent with emission originat- 
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Fig. 7. — Velocity and dispersion of the stellar component derived from the CO bandhead. Contours show contours of the i-C-band image 
with contours separated by a magnitude in surface brightness (as in Fig. 2). Radial Velocities are shown relative to the central velocity of 
-58.9 kms" 1 . 
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Fig. 8. — The LOSVD along the major axis of the innermost com- 
ponent at PA=80°. Panels show the radial velocity, LOS disper- 
sion, h3 and hi components of the LOSVD (from top to bottom). 
Bin centers falling within O'.'l of the major axis are plotted. 
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Fig. 9. — The LOS dispersion values as a function of radius for 
all data points with dispersion errors less than 10 kms -1 . Error 
bars on individual points can be seen in light gray. The thick line 
shows a binned average and standard deviation as a function of 
radius. 

likely the sites of the most recent star formation in the 
nucleus. In §5 we find a total young stellar mass at ages 
< 10 Myr of ~ 10 4 M Q . The central excess has bluer 
optical colors and UV emission but lacks strong Ha or 
Br7 emission, suggesting the young stars there have an 
age of > 10 Myr (Gogarten et al. 2009). 

4. KINEMATICS 

Our NIFS data enable us to derive kinematics for 
both the absorption lines, using the CO bandheads 
from 22900-24000A and the strongest of the H 2 lines at 
21218A. 



ing due to thermal excitation in dense gas as suggested 
by the H 2 2-1 S(l)/l-0 S(l) ratio. 

The distribution of Br7 and Ha visible in Fig. 6 emis- 
sion is primarily concentrated in two clumps, one just to 
the N of the nucleus and one ~1" to the SE. These are 



4.1. Stellar Kinematics 

The velocity and dispersion across the NIFS FOV is 
shown in Fig. 7. At radii of we see a blue-shift 

to the W and red-shift to the E, rotation that is in the 
same direction as the HI rotation seen at larger radii 
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Fig. 10. — Results of kinemetric analysis of the stellar velocity 
field as function of radius. Top - the systemic velocity, middle - 
the position angle, bottom - the amplitude of rotation. Triangles 
in bottom two panels indicate counter-rotation. 



(del Ri'o et al. 2004). At r < 0'.'2, there appears to be 
modest counter-rotation. An alternative view of the ve- 
locity field is seen in Fig. 8 which shows the LOSVD 
profile along a PA=80°. The error bars in the velocity 
range from <1 kms -1 near the center to ~5 kms -1 at 
large radii. The counter-rotation seen in the central ex- 
cess (at r < Of! 2) is clearly significant. However, given 
the nearly face-on orientation of the larger scale of the 
galaxy, this counter-rotation does not necessarily imply 
a large (>90°) difference in the angle of the rotation axes 
of the two stellar components. 

To quantify the position angle (PA) and amplitude of 
the rotation as a function of radius, we used the kineme- 
try program of Krajnovic et al. (2006). This kinemetry 
program uses a generalized version of ellipse fitting to an- 
alyze and quantify the velocity map or other higher order 
moments of the LOSVD. We present the results of the 
kinemetry analysis of the velocity maps in Fig. 10, show- 
ing the systemic velocity, PA, and rotation speed as a 
function of radius. We limited the axial ratio to between 
0.8 and 0.95 (consistent with our ellipse fits to the K 
band image), and fixed the center to the photocenter of 
the dust-corrected K band image. We find an amplitude 
for the central counter-rotation of 3.3 kms -1 (shown in 
triangles in Fig. 10), while the rotation at larger radii 
has a maximum rotation amplitude of 8.4 kms -1 at l'/2. 
The PA of the counter-rotation and rotation at larger 
radii are quite similar, with values of 80-90° (E of N). 
Our observed rotation is at a similar PA to the HI gas, 
which has a rotation axis of 70-90° from 100-300" (del 



Rio et al. 2004). The maximum observed rotation veloc- 
ity of the HI gas is ^45 kms -1 . Recent observations of 
the stellar kinematics along the major axis by Bouchard 
et al. (2010) suggest that increasing counter-rotation rel- 
ative to the HI is seen between radii of 1-20", with ro- 
tation in the sense of the HI again taking over at larger 
radii. Clearly, large-scale IFU kinematics of this galaxy 
would be very interesting. 

A radial profile of the stellar dispersion is shown in 
Fig. 9 for bins with uncertainties in the dispersion of 
<10 kms -1 . While the scatter of values is large at 
most radii, the velocity dispersion clearly declines from 
^40 kms -1 at the center to a mean value of ^33 kms -1 
at 0"6. Our dispersion values are consistent with the 
40±3 kms" 1 determined by Barth et al. (2002) for the 
NGC 404 center using the calcium triplet. Dynamical 
modeling based on the stellar kinematics is detailed in 
§6.1. 

4.2. J?2 Gas Kinematics and Morphology 

The H2 gas kinematics and morphology were derived 
by fitting the 21218A 1-0 S(l) emission line, the strongest 
of the H 2 emission lines in the if-band. In each spaxel, 
we fit this line with a Gaussian to derive a flux, radial 
velocity, and velocity dispersion (after correction for the 
instrumental resolution). The results of these fits arc 
shown in Fig. 11. The physical condition of this gas is 
discussed more in §3.2. We note that the mass of lumi- 
nous excited H 2 gas is very small (< 1 M©) and although 
it likely tracers a larger mass of colder molecular gas, this 
gas is not expected to be dynamically significant in the 
nucleus (Wiklind & Henkel 1990). 

The H2 emission peaks slightly to the SW of the contin- 
uum emission, <~ 0"07 from the centroid of the dust emis- 
sion corrected continuum (see white contours in Fig. 11). 
The FWHM of the distribution is - 0'/4 (6 pc). At 
r < 0"3, it is elongated (q ~ 0.6) with PA~25°. 

The H2 gas velocity field provides a significant contrast 
with the stellar velocity field. It shows a clear rotation 
signature at a position angle nearly orthogonal to the 
stellar rotation. Using the kinemetric analysis program 
(Krajnovic et al. 2006), we find a kinematic PA of ^15° 
out to 0"2, which then twists to PA ~ —10 at larger radii. 
A radial profile of the rotation velocity and dispersion are 
shown in Fig. 12. The rotation peaks at r ~ 0"15 (2.2 pc) 
with an amplitude of 36 kms -1 . Beyond C/3, the rota- 
tion flattens out to an amplitude of 20-25 kms -1 . The 
central dispersion (r = 0"05) is 39 kms -1 , which likely 
includes unresolved rotation. At r = O'.'l it drops to 
33 kms -1 , then slowly declines to ^25 kms -1 at larger 
radii. Dynamical modeling of the gas kinematics are pre- 
sented in §6.2. 

The axial ratio of the intensity distribution and our 
gas dynamical model suggest an inclination of ~30° at 
r < 0"3, with smaller inclinations at larger radii. As- 
suming this geometry, the H 2 disk appears to be rota- 
tion dominated at all radii with V TOi /a <~ 2 near the 
center. H 2 emission line disks are commonly observed in 
nearby AGN (Neumayer et al. 2007; Riffel et al. 2008; 
Hicks et al. 2009), although the size of these disks are 
typically much larger than in NGC 404. For instance, in 
the Hicks et al. (2009) sample of 11 nearby AGN, they 
find H 2 disks traced by the 21218A line with half-width 



The NGC 404 Nucleus 



11 



X Offset [pc] 
40 10 



1 




-20 



X Offset [pc] 
-10 10 



50 u 

t/3 



30 



10 | 
o 

-10 | 

| 

-30 | 

as 

-50 * 





1.0 




0.5 






0.0 


Q 






-0.5 
-1.0 I 




55 



44 8 



33 eM, 

pi 
o 

22 'I 
u 

1/3 

ii s 
o 



-1.0 -0.5 0.0 0.5 
X offset ["] 



1.0 



-1.0 -0.5 0.0 0.5 
X offset ["] 



1.0 



Fig. 11. — Velocity and dispersion of the molecular hydrogen data determined from the 1-0 S(l) line at 21218A. Black contours indicate 
the surface brightness of the A'-band image with contours separated by a magnitude in surface brightness. White contours show the H2 
total intensity 5, 10, 20, 30, 40, 50, 70, and 90% of the peak. Radial velocities are shown relative to the central stellar radial velocity (-58.9 
kms -1 ). 
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Fig. 12. — The rotation velocity (black line) and dispersion (blue 
dashed line) derived from H2 21218A line using the kinemetry pro- 
gram of Krajnovic et al. (2006). 



half-maximum values ranging from 7-77 pc (vs. 3 pc in 
NGC 404). Most of these disks have V TO tJa values of —1, 
although a couple have V TOt /a > 2 at larger radii. 

5. STELLAR POPULATIONS 

In this section we use the optical spectrum of the 
NGC 404 nuclear star cluster to constrain its formation 
history. We first use the optical emission lines to con- 
strain the reddening and derive a maximum current SFR 
of 1.0 x 10~ 3 M0/yr. Fits to the absorption line spectrum 
show that the nucleus has a range of stellar ages, with 
stars ~1 Gyr in age contributing more than half the light 
of the spectrum. Next, we consider the range in mass- 
to- light ratios (M/ Ls) allowed by the stellar populations 
fits. Finally, we show that the bulge has a distinctly older 
stellar population than the nucleus. 



5.1. Optical Emission Line Information 

Before discussing our fits to the optical absorption line 
spectrum, we first use the optical emission lines to esti- 
mate the reddening, current metallicity and star forma- 
tion rate (SFR). We note again that our line fluxes and 
ratios agree very well with those in Ho et al. (1997). 

After correcting for the underlying stellar population 
using our population fits described below, we find a 
Balmer decrement (F(Ha)/F(Hj3)) of 3.59. For an AGN 
the expected Balmer decrement is 3.1 while for an HII re- 
gion it 2.85 (Osterbrock 1989). Given the ambiguous na- 
ture of the nuclear emission, this suggests extinctions of 
A v = 0.43-0.73. As noted in §3.1, the F547W-F8UW 
color-map indicates extinctions of Ay < 2 covering about 
half the nucleus, thus these Balmer decrement extinction 
values seem reasonable. We note that the patchy distri- 
bution of the Ha emission (see Fig. 6) is significantly dif- 
ferent from the stellar light distribution, thus the Balmer 
decrement extinction may not exactly match the extinc- 
tion of the stellar light. 

If we assume that the emission is primarily powered 
by star-formation, then we can estimate the metallic- 
ity based on strong-line indicators. Using the metallic- 
ity calibration of the [Oil] /[Nil] ratio by Kewley & Do- 
pita (2002), we find a metallicity of 12+log(0/H)=8.84, 
roughly solar. We can also estimate the gas-phase 
metallicity from the Tremonti et al. (2004) metallicity- 
luminosity relationship; assuming Mb ~ —16.75 (using 
values from Tikhonov et al. 2003), we get a metallic- 
ity 12+log(0/H)=8.33. These results suggest a solar or 
slightly sub-solar current gas-phase metallicity. 

Finally, we can place an upper limit on the current 
star-formation rate in the nucleus using the Ha lumi- 
nosity. Correcting for extinction, we get an Ha flux 
of 10.8xl0~ 14 ergss~ 1 cm -2 corresponding to Ln a = 
1.2 x 10 38 ergss" 1 . This gives an upper limit to the cur- 
rent star formation rate of 1.0 x 10 -3 M Q yr _1 using the 
Kennicutt (1998) relation. The value is an upper limit 
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as some of the Ha flux may originate from AGN activity. 
The radio continuum detection of 3 mJy at 1.4 cm by 
(del Rio et al. 2004) implies a SFR of -2x 10" 4 M© yr" 1 
if it is the result of recent star formation. 

The UV spectra of Maoz et al. (1998) in an 0"86 aper- 
ture also suggests the presence of O stars in the nucleus. 
We correct their estimate of the number of O stars to a 
distance of 3 Mpc and for extinction of Ay = 0.5 — 1.0. 
This gives a range of 17-170 O stars, suggesting a to- 
tal stellar mass of 4000-40000 M Q assuming a Kroupa 
et al. (1993) IMF. The low end of these estimates are 
consistent with the 10 3-4 M© of recent stars implied by 
the Ha and radio SFRs assuming formation over the last 
-10 Myr. 

5.2. Nuclear Cluster Stellar Populations 

In this section we fit the integrated nuclear optical 
spectrum to determine the stellar populations of the 
NSC. We find that stars with a range of ages are present, 
but that stars 1-3 Gyr make up >50% (6 x 10 5 M Q ) of 
the NSC mass. This dominant population has a solar or 
slightly subsolar metallicity ([Z]=-0.4). About 10 4 M 
of young stars are also present. These results are con- 
sistent with other studies of the NGC 404 nucleus (Cid 
Fernandes et al. 2004; Bouchard et al. 2010). 

We fit the flux-calibrated spectrum between 3740- 
6200A, avoiding the second order UV contamination at 
wavelengths above 6200A. For stellar population mod- 
els we use both (1) the Bruzual & Chariot (2003) mod- 
els (BC03) which are based on Padova stellar models 
(Bertelli et al. 1994; Girardi et al. 2000) and the STELIB 
library of stellar spectra (Le Borgne et al. 2003) and (2) 
the updated Chariot & Bruzual 2007 (CB07) version of 
these models incorporating the AGB models of Marigo 
& Girardi (2007) and Marigo et al. (2008) obtained from 
Stephane Chariot. These models provide spectra for sin- 
gle age simple stellar populations (SSPs) over a finely 
spaced age grid and at metallicities of Z=0.0001, 0.0004, 
0.004, 0.008, 0.02, and 0.05 at a spectral resolution simi- 
lar to our observations. We find that typically the CB07 
models more accurately describe the data; they differ 
from the BC03 models by up to 10% across the wave- 
length range we are considering here. 

To remove known problems in the wavelength calibra- 
tion of the STELIB library (Koleva et al. 2008) and 
any wavelength issues in our own spectrum, we deter- 
mined the relative velocity between our spectrum and 
the SSP models every 100A. Variations of up to ±35 
kms -1 around the systemic velocity were found and cor- 
rected. We also made a flux calibration correction of 
< 3% across the spectrum that was determined from 
residuals of our fits to the nuclear clusters in NGC 404, 
NGC 205, M33, and NGC 2403 - this is described in 
more detail below. We fit the models to our spectrum us- 
ing Christy Tremonti's simplefit program (Tremonti et al. 
2004), which uses a Levenberg-Marquardt algorithm to 
perform a \ 2 minimization to find the best-fit extinction 
(using the prescription of Chariot & Fall 2000) and scal- 
ings of the input model spectra. The code excludes the 
regions around expected emission lines from the fit. 

We started our fitting by considering just single 
age stellar populations. Over the full range of 3740- 
6200A, the best-fitting SSP is the CB07 model with 



Agc=1.14 Gyr, Z=0.004, and an Ay = 0.85 mags. The 
reduced \ 2 1S 21.8 corresponding to typical residuals 
of 3%. This fit is shown in Fig. 13 as the blue line. 
The residuals are dominated by a mismatch between the 
large-scale shape of the nuclear cluster continuum with 
the SSP model. To mitigate these effects and get a bet- 
ter sense of the metallicity, we also tried fitting just the 
metal line rich region from 5050-5700A. The best fit is 
again a CB07 model with an age of 1.14 Gyr, but with 
Z=0.008 and A v = 0.91. This metallicity and redden- 
ing are quite similar to the values expected based on the 
emission lines and the mass-mctallicity relationship. 

Previous studies have found significant evidence for 
multiple stellar populations in nuclear star clusters (Long 
et al. 2002; Rossa et al. 2006; Walcher et al. 2006; Seth 
et al. 2006; Siegel et al. 2007). While most of these stud- 
ies have focused on spiral galaxies, there is also direct 
evidence that the NGC 404 nuclear star cluster contains 
multiple stellar populations; Maoz et al. (1998) finds ev- 
idence for massive O stars, suggesting a population with 
a g e ^10 Myr, while based on the color of the nucleus 
and its spectrum, older stars are clearly present as well. 
Also, as discussed in §3.1, we expect contamination of 
—30% from the bulge in our MMT spectrum. 

We tried a number of approaches to our multi-age fits. 
All of these involved fitting a small subsample (—10) of 
the SSPs with ages ranging from 1 Myr to 20 Gyr (sec 
Walcher et al. 2006; Koleva et al. 2008, for a more com- 
plete discussion of stellar population synthesis) . We tried 
three types of models with both the CB07 and BC03 
spectra: 

1. Single-metallicity models using the 10 default ages 
in the simplefit code (5, 25, 101, 286, 640, 904, 
1434, 2500, 5000, 10000 Myr) or the 14 age bins 
used by Walcher et al. (2006) (1, 3, 6, 10, 30, 
57, 10i, 286, 570, 1015, 3000, 6000, 10000, 20000 
Myrs). 

2. Inspired by observations of the Sgr dSph nucleus 
(a.k.a. M54; Monaco et al. 2005; Siegel et al. 2007) 
we also tried fitting models with chemical evolution 
from low metallicity (Z=0.0004 / [M/H]=-1.7) at 
old ages (13 Gyr) to solar metallicity for young ages 
(following Fig. 2 from Siegel et al. 2007). Specifi- 
cally we used a model with Z=0.02 at ages of 1, 10, 
50, 101, 286, 570, 1015, and 2500 Myr plus Z=0.008 
models at 5000 and 10000 Myr and Z=0.0004 mod- 
els at 13000 Myr. 

3. Assuming a continuous or exponential star forma- 
tion history with chemical evolution similar to the 
previous item. 

Of all the single metallicity models (model 1), the best 
fit is obtained with the CB07 spectra using the default 
ages with Z=0.008 yielding a reduced \ 2 of 4.3 and typi- 
cal residuals of 1%. This is a vast improvement over the 
best fit SSP model. The fit is dominated by interme- 
diate age stars, with 72% of the light in the 1.434 Gyr 
SSP. Many other models also produce nearly as good fits 
as this one including models with solar and super-solar 
metallicity. 

All the best-fitting {x 2 re d < 5) models have fairly sim- 
ilar age distributions with 50-85% of the light coming 
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Fig. 13. — (Top panel -) Normalized spectrum of the NGC 404 nucleus (black) with the best-fitting single age (blue) and multi-age (red) 
fits to the spectrum. Vertical lines show emission line regions that were excluded from the fit. Gray spectra indicate the different aged 
SSPs contributing to the multi-age fit. (Bottom panel) - residuals of the best fit single age (blue) and multi-age (red) fits. Also shown is 
the flux calibration correction described in the text. 



TABLE 3 

Best-Fit NSC Stellar Population Model 



Age [Myrj 


1 


10 


50 


101 


286 


570 


1015 


2500 


5000 


10000 


13000 


Z 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.008 


0.008 


0.0004 


Light Fraction 


0.013 


0.090 


0.000 


0.000 


0.157 


0.000 


0.277 


0.239 


0.166 


0.058 


0.000 


Mass Fraction 


0.000 


0.003 


0.000 


0.000 


0.035 


0.000 


0.152 


0.317 


0.316 


0.177 


0.000 



from ages between 600 and 3000 Myrs, 5-6% of the light 
coming from ages <10 Myr, and 0-25% of the light com- 
ing from ages >5 Gyr. The Ay values range from 0.84 to 
1.16 and the implied V — I colors range from 0.84-1.04. 

The spectrum is not well fit by a constant star forma- 
tion history (model 3), giving a reduced x 2 of 22. The 
constant SFH clearly has too many young stars, with the 
continuum mismatched in a way that cannot be com- 
pensated for by the reddening. We note that a constant 
SFH does provide good fits to the spectra of nuclei in 
many late- type galaxies (Seth et al. 2006; Walcher et al. 
2006). An exponentially declining SFH with a timescale 
of ^8 Gyr does a better job of matching the overall con- 
tinuum of the observations, but still yields a reduced x 2 
double that of our best fitting models. 

The best-fit model is obtained from model 2, which 
contains metallicity evolution. The reduced x 2 of the 



best fit is 4.23, slightly better than any of the single 
metallicity models. This model is shown in red in Fig. 13 
with the individual SSP components shown in gray in the 
top panel. The best-fit model has significant contribu- 
tions from many of the individual ages, these are shown 
in Table 3. It is once again dominated by 1-3 Gyr old 
stars, which make up half the mass of the NSC. This 
appears to be a robust result, with all our best-fit mod- 
els containing a similar result. The luminosity weighted 
mean age is 2.3 Gyr, the mass weighted mean age is 
4.3 Gyr. The best fit Ay is 0.88 and the M/Ls are 0.89, 
0.73, and 0.34 in the V-, I-, and #-bands. NGC 404 
does not seem to have a dominant old metal-poor popu- 
lation like that seen in M54 (Monaco et al. 2005). From 
the M/L analysis presented in §5.3, we find that the 
maximum contribution for an old (13 Gyr) metal-poor 
(Z=0.0004) population that keeps the reduced x 2 below 
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5 is ~17% of the luminosity in the U-band and ^35% of 
the total mass. 

We can compare our best-fit model to other observed 
properties of the nucleus. The integrated colors of the 
best-fit model are I - #=1.59 and V - 7=0.96. An es- 
timate of the I — H <~ 1.7 (Table 2) comes from the SB 
profiles, while the F547W - F8UW ~ V - I ~ 0.95 
comes from the unreddened W portion of the nucleus in 
Fig. 2. Both of these quantities need to be corrected for 
foreground extinction of E(I — H) and E(V — I) of 0.06 
and 0.08 magnitudes; after this correction the measured 
integrated colors agree with the best-fit spectral synthe- 
sis model to within 10%. 

The best-fit model also suggests that the young stars 
(< 10 Myr) make up 6%/0.2% of the light/mass in H- 
band and 8%/ 0.3% of the light /mass /-band. This gives 
a mass of of about 1.1 x 10 4 M Q in young stars, in very 
good agreement with Ha SFR integrated over 10 Myr 
and the estimates from the UV spectra (see discussion in 
§5.1; Maoz et al. 1998). 

Finally, we compared the spectrum of NGC 404 to 
other nuclear star clusters for which we took spectra on 
the same night, including M33, NGC 205, NGC 2403, 
and NGC 2976. Both M33 and NGC 205 do better 
than any single stellar population model at matching the 
NGC 404 spectrum, with reduced \ 2 values of 9.3 and 
15.0. The stellar populations in the M33 nucleus have 
been previously discussed by Long et al. (2002). They 
find the spectrum requires a large contribution of stars 
with an age of ~1 Gyr and a lesser contribution of stars 
at 40 Myr. This dominant 1 Gyr population is consis- 
tent with our findings above. In our initial fits of these 
nuclear star clusters we noted that the residual pattern 
was quite similar regardless of the details of the fit, with 
the most dramatic feature being a change of a few per- 
cent between 4700-4800A. Our flux calibration errors are 
expected to be at this level based on the residuals from 
the flux calibrators. Given the similarities of the residu- 
als in all of the nuclear spectra we medianed the best fit 
model 1 residual spectra from NGC 205, M33, NGC 404, 
NGC 2403, then smoothed it and used it as a correction 
to our flux calibration. The correction is shown as the 
green line in the residual plot in Fig. 13. This change im- 
proved the x 2 01 the single-stellar population and multi- 
age fits substantially at all metallicities but had little 
effect on the best-fit parameters. 

The stellar population of NGC 404 has been previously 
studied by Cid Fernandes et al. (2004) who analyzed the 
blue spectrum (3400-5500A) of the NGC 404 nucleus as 
part of a study of a large number of LINER nuclei. Using 
empirical templates of five nuclei ranging from young to 
old they find a best fit for NGC 404 of 78.5% of light 
in an intermediate age template (for which they used a 
spectrum of NGC 205) and 21.5% from old templates, 
with an Ay of 0.90 and residuals of 6.6%; this residual is 
much larger than for our fits (with <~1% residuals) due 
primarily to the lower S/N of their spectrum. Their fit 
agrees quite well with our finding of a ~1 Gyr old pop- 
ulation in the nucleus and our best-fit reddening. Un- 
like us, they find no evidence for young stars, probably 
because the NGC 205 template spectrum used for in- 
termediate ages also includes a young stellar component 
(e.g. Monaco et al. 2009). The nuclear stellar popula- 



tions have also been recently studied by Bouchard et al. 
(2010) using a different fitting code and stellar library. 
Their best-fit four-component model matches ours very 
well, with 4% of the light in stars younger than 150 Myr, 
20% at 430 Myr, 62% at 1.7 Gyr, and 14% at 12 Gyr. 

5.3. Nuclear Cluster Mass to Light Ratio 

The multi-age population fits to the NGC 404 nuclear 
spectrum are degenerate, with many possible combina- 
tions of ages and extinctions giving similar quality fits. 
This translates into an uncertainty in the stellar popu- 
lation M/L, with the presence of more young stars de- 
creasing the M/L and the presence of more old stars in- 
creasing it. For comparison with the dynamical models 
(§6.1), we'd like to quantify the range of possible M/Ls 
allowed by spectral synthesis fits. 

To try to quantify the possible variation in the M/L, 
we took the best-fit model described in the previous sec- 
tion (model 2, with CB07 and current solar metallicity) 
and fixed the age of each individual SSP (e.g. just the 
10 Gyr population) to have between 0-100% of the flux 
of the unconstrained best-fit model, while the extinction 
and contribution of the rest of the SSPs were allowed to 
vary. We then took all 'acceptable' fits (defined below) 
and examined the range in M/L. Because we have SB 
profiles in both I- and Ti-bands (Fig. 3), we focus on the 
M/Ls in those bands. 

We choose to consider all fits with reduced \ 2 values 
of < 5 as acceptable, a value 18% higher than our best 
fitting models. This choice is informed by the \ 2 proba- 
bility distribution; our fits have ~660 degrees of freedom 
and ignoring any model errors, the 99% confidence limit 
for this distribtion is a A\ 2 of 13%. Unfortunately, since 
our best fit has a reduced \ 2 is 4-2, there must be signif- 
icant deficiencies in the models or flux calibration at the 
1% level. Nonetheless, the range of M/Ls with reduced 
X 2 < 5 should reflect the a conservative estimate of the 
acceptable values, roughly equivalent to 3— a error bars. 

We find that the range of acceptable M/Ls is 0.68- 
1.22 in I band and 0.29-0.61 in H band. We note that 
the best-fit M/L of 0.73 in I band and 0.34 in H band 
is near to the lower range of these values, suggesting a 
larger degeneracy in the mass of old stars than for the 
younger populations. Our best-fit values are consistent 
with with the dynamical M/Lj = 0.70 ± 0.04 derived in 
§6.1. 

5.4. Bulge Stellar Populations 

In Cid Fernandes et al. (2005), they analyze the radial 
variations in the NGC 404 spectrum and find increas- 
ing age with radius, with the old and intermediate age 
templates having equal weight at radii between 5 and 10 
arcseconds. We use our MMT spectra of NGC 404 to 
extract a spectra at radii between 5" and 10 "from the 
nucleus. From our profile fits (Fig. 3), we expect that the 
bulge component should dominate this spectrum. This 
spectrum has a peak S/N of 70 around 5000A dropping 
to 25 at 3700A. 

The models that best fit the NGC 404 bulge spec- 
trum are the same sets of models as the nuclear spec- 
trum (i.e. model 1 with Z=0.008 and model 2). How- 
ever, like Cid Fernandes et al. (2005) we find that the 
bulge has a significantly older population with all the 
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best fits having >50% of the light in an old (> 5 Gyr) 
component. The best fits have Ay between 0.16 and 0.24 
which corresponds well to the Ay of 0.19 expected for the 
foreground extinction (Schlegel et al. 1998) and is consis- 
tent with the dust absorption being primarily restricted 
to the inner 5"of the galaxy as suggested by HST color 
images. The best fits all have residuals of ^2.6%, corre- 
sponding to reduced % 2 of 2.0. All the fits have ~2% of 
the light in a young (<100 Myr) component and a signif- 
icant contribtution (~30%) from a ~1 Gyr component. 
This analysis clearly shows a significant population dif- 
ference between the bulge and nuclear spectrum which 
supports the separation of the SB profiles into separate 
components in §3. 

As all the models have nearly identical fits, we infer 
M/L ratio using model 2, which best fit the nuclear star 
cluster spectrum. We find a M/L of 1.69, 1.28, and 0.63 
in the V-, I-, and //-bands. 

Two contemporaneous papers suggest that a dominant 
old population extends out to larger radii as well. At r ~ 
25", Bouchard et al. (2010) find two-thirds of the mass is 
an old component using spectral synthesis. And using a 
deep color-magnitude diagram at a radius of ^160"(near 
the disk-bulge transition), (Williams et al. 2010) have 
found ^80% of the stellar mass is >10 Gyr old. The 
nucleus, with its dominant intermediate age population, 
is therefore quite distinct from the rest of the galaxy. 

6. DYNAMICAL MODELING 

In this section we present dynamical modeling of the 
central mass distribution of NGC 404 using both the 
stellar (absorption line) and gas (H2 emission) kinemat- 
ics. From modeling of the stellar kinematics, we find a 
M/Lj = 0.70 ± 0.04. We can use this to derive a total 
mass for the nuclear star cluster by combining this esti- 
mate with the total luminosity of the NSC and central 
excess (Table 2) to get a mass of 1.1 ±0.2 x 10 7 M Q . The 
error on the NSC mass includes additional uncertainty 
due to the possible variation in M/L of the central ex- 
cess. We note that the dynamical M/Lj matches the 
best-fit M/Li derived from the stellar population fits to 
within 10%. 

We can also constrain the presence of a MBH using 
the dynamical models. The two methods give inconsis- 
tent results for the BH mass, with the stellar kinematics 
suggesting an upper limit of < 1 x 10 5 M©, while the gas 
models are best fit with a MBH of mass 4.5±2' 5 x 10 5 M Q 
(3<T errors). These measurements assume a constant 
M/L within the NSC; taking into account possible vari- 
ations of M/L at very small radii, the BH mass esti- 
mates/upper limits can increase by up to 3 x 10 5 M Q . If 
there is a MBH in NGC 404, its mass is constrained to 
be < 10 6 M Q , smaller than any previous dynamical BH 
mass measurement in a present-day galaxy. 

6.1. Stellar Dynamical Model 

We constructed a dynamical model to estimate the 
mass and M/L of the nuclear star cluster and of a pos- 
sible central supermassive black hole inside it. The first 
step in this process is developing a model for the light 
distribution. To parametrize the surface brightness dis- 
tribution of NGC 404 and deproject the surface bright- 
ness into three dimensions, we adopted a Multi-Gaussian 
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Fig. 14. — x contours describing the agreement between the 
JAM dynamical models and the NIFS integral-field observations 
of Vrms- The models optimize the three parameters Mbh, 0z and 
M/L but the contours are marginalized over the M/L, which is 
overplotted as dashed contours with labels. Dots show the grid of 
models. The three lowest \ 2 contour levels represent the Ax 2 = 1, 
4 and 9 (thick red line) corresponding to 1, 2, and 3cr confidence 
levels for one parameter. 




1.0 

R (arcsec) 

Fig. 15. — The observed NIFS profiles of Vrms, biweight averaged 
over circular annuli (data points), are compared to JAM model pre- 
dictions of the same quantity, for different Mbh = (lowest solid 
line), 1, 3, 5, and 10x10 s Mq (top dashed line). All models have 
M/Lj = 0.70 Mq/Lq and /3 Z = 0.5 as the best fitting solution. 

Expansion (MCE; Emsellem et al. 1994). The MCE fit 
was performed with the method and software 10 of Cap- 
pellari (2002). We fit a 2-D form of the J-band data as 
presented in §3: at r < 5", we use the WFPC2 /-band 
image with reddened regions masked out, while outside 
of 5", we use the ACS data (which is saturated in the nu- 
cleus) . Although the use of NIR data would be preferable 
given the presence of dust in the nucleus, the uncertainty 
in the NIFS PSF and the poor resolution and sampling of 
the available HST NICMOS data makes the F814W data 
the best available image. The MCE fit takes into account 
the PSF: we use a Tiny Tim PSF (Krist 1995). In the fit 
we forced the observed axial ratios q' of the MGE Gaus- 
sians to be as big as possible, not to artificially constrain 
the allowed inclination of the model. The PSF-convolved 
MGE represents a reasonable match to the F814VF-band 
images with some mismatch seen in the inner parts due to 
dust. We note that modeling results based on NICMOS 
F160W mass models are consistent with those presented 
here. 

10 Available from http://purl.org/cappellari/idl/ 
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The primary uncertainty in the determination of Mbh 
results from uncertainties in the central stellar mass pro- 
file. The F8UW luminosity at r < 0705 (0.7 pc) is 
~ 1.3 x 10 6 L Q , while the MGE model has a luminos- 
ity of ~ 4 x 10 5 L Q within a sphere of the same radius. 
Translating this luminosity into a stellar mass density is 
complicated given the possible presence of dust, stellar 
population gradients, non-thermal emission and possible 
nuclear variability. Furthermore, the available kinemat- 
ics sample only the very nuclear R < l'/5 region of the 
galaxy. This spatial coverage is not sufficient to uniquely 
constrain the orbital distribution and Mbh in NGC 404 
(Shapiro et al. 2006; Krajnovic et al. 2009). This complex 
situation does not justify the use of general brute-force 
orbital superposition dynamical modeling methods (e.g. 
Richstone & Tremaine 1988; van der Marel et al. 1998). 
This implies we will not be able to uniquely prove the 
existence of a BH in this galaxy, but we can still explore 
the ranges of allowed Mbh by making some simplifying 
observationally-motivated assumptions. A similar, but 
less general approach was used in other studies of BHs 
in NSCs from stellar kinematics (Filippenko & Ho 2003; 
Barth et al. 2009). We note that because we kinemat- 
ically resolve the NSC and can determine its M/L, we 
can provide stricter BH upper limits than those provided 
by these papers, in which the NSC kinematics were un- 
resolved. 

The adopted dynamical model uses the Jeans 
Anisotropic MGE (JAM) software 10 which implements 
an axisymmetric solution of the Jeans equations which 
allow for orbital anisotropy (Cappellari 2008). Once a 
cylindrical orientation for the velocity ellipsoid has been 
assumed, parameterized by the anisotropy parameter 
(3 Z = 1 — /a ft, the model gives a unique prediction for 

the observed second velocity moments V lms = y/V 2 + <r 2 , 
where V is the mean stellar velocity and a is the veloc- 
ity dispersion. The JAM models appear to provide good 
descriptions of integral-field observations of large sam- 
ples of real lenticular galaxies (Cappellari 2008; Scott 
et al. 2009) and give similar estimates for BH masses as 
Schwarzschild models (Cappellari et al. 2010). 

The luminous matter likely dominates in the extreme 
high-density nucleus of NGC 404. For this dynamical 
model we assume light traces mass with a constant M/L. 
For the model we assumed a nearly face-on inclination of 
i = 20° (i = 90° being edge on) as indicated by the regu- 
lar HI disk kinematics and geometry (del Rio et al. 2004), 
and further supported by the nearly circular isophotes at 
small radii. However, unlike the gas dynamical model, 
our results do not depend strongly on this inclination 
assumption. The model has three free parameters: (i) 
The anisotropy /3 Z , (ii) the mass of a central supermas- 
sive black hole Mbh and (iii) the /-band total dynam- 
ical M/L. To find the best fitting model parameters 
we constructed a grid for the two nonlinear parameters 
(P z , Mbh) values and for each pair we linearly scaled the 
M/L to match the data in a \ 2 sense. The resulting 
contours of x 2 , marginalized over M/L are presented in 
Fig. 14. They show that with the assumed stellar mass 
distribution, the models suggest a Mbh < 1 x 10 5 M 
with (3 Z ~ 0.5 and M/L = 0.70 ±0.04 M /L Q in /-band 
(3cr levels or Ax 2 =9). A corresponding data-model 
comparison is shown in Fig. 15. The dynamical M/Lj 



(which includes a correction for foreground extinction) is 
quite close to our best-fit stellar population estimate of 
M/Li = 0.73. 

As noted above, the largest uncertainty in our MBH 
measurements comes from uncertainty in the central stel- 
lar mass density hidden by the assumption of a constant 
M/L. The bluer color in the F547M- F8UW color im- 
age (Fig. 2) and the presence of UV emission at r < 0'.'05 
suggests a possible decrease in M/L which could mask 
the presence of a possible BH. Because our optical spec- 
troscopy does not resolve the central stuctures, we can 
only use the HST imaging data to try constraining possi- 
ble variations in the M/L. We now determine a limit on 
how much the central M/L may vary by quantitatively 
examining the color difference between the NSC and the 
central pixels. We measured the NSC colors just to the 
W of the nucleus and find a F330W - F547M (U - V) 
of 0.48 and an .F547 - F814VF (V - I) of 0.87 after cor- 
rection for foreground reddening. The colors at r < 0'.'05 
are both about 0.1 magnitudes bluer than the NSC col- 
ors. In V — I the presence of a ~10 Myr population will 
give the largest reduction in M/L with the least affect 
on the color, since these populations are dominated by 
red supergiants at long wavelengths (e.g. see Fig. 4 in 
Seth et al. 2008b). For a solar metallicity population, the 
V — I color at 10 Myr is about 0.5 (Girardi et al. 2000), 
meaning the observed nuclear color difference is consis- 
tent with about 30% of the light in young stars. This 
would result in a ^30% reduction in the M/L, as the 
M/L for a 10 Myr population is very small. However, 
this large fraction of a young population would imply a 
larger change in U — V color than observed, so 30% rep- 
resents an upper limit on the reduction in M/L. If we 
take the /-band luminosity at r < 0'.'05 of 1.3 x 10 6 L Q 
with our default M/L of 0.70, we get a mass at r < 0'.'05 
of 9 x 10 5 M . Reducing the M/L by 30% at the center 
relative to the NSC would therefore reduce the central 
stellar mass and increase the BH mass by <~ 3 x 10 5 M . 
Therefore, even including the effects of a possible decline 
in the central M/L, the BH mass is constrained to be 
< 4 x 10 5 Mq by the stellar dynamical model. Future 
HST spectroscopy that can resolve the stellar popula- 
tions gradient in the nucleus should produce a more ac- 
curate mass model that will enable us to robustly test 
the presence of a MBH in NGC 404. 

6.2. Gas Dynamical Model 

In addition to the stellar dynamical model we con- 
structed a dynamical model to reproduce the kinematics 
of molecular hydrogen emission (§4.2) to obtain an ad- 
ditional constraint on the presence of a possible central 
BH. We follow the modeling approach of Neumayer et al. 
(2007) that assumes the gas to be settled in a thin disk 
in the joint potential of the stars and a putative BH. We 
assume the gas is moving on circular orbits and we do not 
include a pressure term in the gas model, as the velocity 
dispersion of the gas is small. We have tested the effect of 
an additional pressure term and find it to be very small, 
well below our derived errors. The stellar potential is 
fixed using the MGE parameters and the stellar M/Lj 
derived §5.2, and thus suffers from the same uncertainty 
in the central stellar mass profile. 

The gas velocity map appears quite complex with sud- 
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Fig. 16. — Velocity field of the best-fitting H2 dynamical model (middle, top), with a black hole mass of Mbh = 4.5 X 10 5 Mq and a 
disk inclination of 37°, in comparison to the symmetrized data (top left). The velocity residual (data-model) is shown in the right panel. 
The velocity curves in the bottom panels are extracted along the line of nodes (overplotted to the velocity maps), and represent the peak 
velocity curves. The diamonds correspond to the model velocity curve, while the crosses correspond to the data. 

the stellar model and consistent with inclinations derived 
by del Rio et al. (2004) from HI. The dynamical model 
has two free parameters: (i) The inclination angle ol the 
inner (r< 0'.'2) H2 disk i and (ii) the mass of a central 
supermassive black hole A/eh- To find the best fitting 
model parameters we constructed a grid for the two pa- 
rameters (i,Mbh) and matched the data (velocity and 
velocity dispersion) minimising \ 2 . Figure 17 shows the 
contours of x 2 an d indicates the best fit model param- 
eters i = 37° ± 10° and M BH = 4.5±l;g x 10 5 M (3-cr 
errors). The best-fitting model for the H2 velocity field 
is compared to the data in Fig. 16. We note that the 
best-fit inclination is consistent with the axial ratios of 
the emission (§4.2). 

As Figure 17 shows, the inclination angle i and the cen- 
tral mass Mbh are strongly coupled, since the amplitude 
of the rotation curve is proportional to \/Mbh X sin(i). 
This degeneracy in the gas model leads to the relatively 
large uncertainties in inclination angle and black hole 
mass. However, the gas kinematics cannot be repro- 
duced without a central black hole, even for an almost 
edge-on central inclination angle. The best-fitting black 
hole mass is significantly different from the non-detection 
of < 1 x 10 5 Mq derived from the stellar kinematics. 
One possible source of this discrepancy would be non- 
gravitational motions that are not accounted for by the 
gas-kinematical model. Another possible problem is the 
assumption of a thin gas disk, although thicker gas dis- 
tributions provided worse fits to the gas kinematic data. 
We note that both determinations use the same MGE 
mass model, and thus the uncertainty due to the central 
M/L applies to these results as well. This could increase 
the BH mass estimates by as much as 3 x 10 5 Mq. 



4 6 

M BH (10 5 M Q ) 

Fig. 17. — Constraining the mass of the central BH: the figure 
indicates the grid of models (in black hole mass, Mbh, and disk 
inclination) that was calculated, and the contours show x 2 tn the 
vicinity of the best-fit dynamical models for matching the H2 kine- 
matics. The minimum x 2 model is at Mbh ~ 4.5 X 1O 5 M0 and 
a disk inclination of 37°. The contours indicate the 1, 2, and 3 a 
confidence levels, respectively (see text for details). 



den changes in position angle of the major kinematic 
axis and some of the kinematic features might be due 
to non-gravitational motions driven by the LINER nu- 
cleus. However, as described in §4.2, there also appears 
to be clear rotation in the central 0'.'4 x 0'.'4 with a major 
rotation axis of ~ 15°, and a twist to slightly negative 
PAs at larger radii. We interpret this twist in the line-of- 
nodes as a change in inclination angle of the molecular 
gas disk and model the kinematics using a tilted ring 
model. The inclination angle of the gas disk for radii 
> 0'.'2 is fixed to the value of 20° that is also used for 
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Fig. 18. — Scaling relations for NSCs and MBHs in early-type 
galaxies. Data for MBHs is from Haring & Rix (2004), while for 
the NSCs it is from Cote et al. (2006) with mass estimates from 
Seth et al. (2008a). The solid line shows the joint MBH/NSC fit 
from Ferrarese et al. (2006), while the dashed line shows the fit to 
the MBH data from Haring & Rix (2004). Just the bulge mass is 
plotted for NGC 404, the total galaxy mass would be ~30% larger. 



7. DISCUSSION 

7.1. Nuclear Star Cluster: Scaling Relationships and 
Formation 

We compare the NGC 404 NSC to the scaling rela- 
tionships between galaxy/bulge mass and BH/NSC mass 
in Fig. 18. We have derived the mass of the NSC us- 
ing the dynamical Mj L determined from the NIFS stel- 
lar kinematics and the luminosity of the NSC (and cen- 
tral excess) above the underlying bulge. The NGC 404 
NSC mass is a large fraction (1-2%) of the bulge/galaxy 
mass, making it fall significantly above the relationship 
for MBHs and early-type Virgo NSCs given by Ferrarese 
et al. (2006), however, there are several other similar out- 
liers in the Virgo NSC sample. Using the Ferrarese et al. 
(2006) NSC mass vs. a relation also gives a mass nearly 
an order of magnitude smaller than the measured NSC 
mass. 

The NSCs high mass relative to its galaxy raises the 
question of whether our interpretation of this component 
as an NSC is correct. As an alternative, we could be 
seeing some form of "extra-light component" (Hopkins 
et al. 2009; Kormendy et al. 2009) or an "extended nu- 
clear source" (Balcells et al. 2007) that are typically seen 
on somewhat larger scales than NSCs in coreless ellipti- 
cal, lenticular, and early-type spirals. The mass fraction 
of —1% places the NGC 404 NSC in the regime of overlap 
between NSCs and the extra-light components (Fig. 46 
of Hopkins et al. 2009). 

Using data collected in Seth et al. (2008a), we plot the 
size and mass of clusters in Fig. 19, and find that the NSC 
in NGC 404 is typical of those in other early- type galax- 
ies. The NGC 404 NSC is also typical in the mass-density 
plane, and the NSC Sersic index of ^2 is similar to that 
found for the M32 and Milky Way NSC in (Graham & 
Spitler 2009). The morphology of the NGC 404 NSC 
is therefore completely consistent with other NSCs. Its 
mass and density is also similar to the lowest mass extra- 
light components (Hopkins et al. 2009, Fig. 45), and fol- 
lows the size-surface brightness scaling of the extended 
nuclear components (Fig. 3 of Balcells et al. 2007). In 
Hopkins et al. (2009), they argue that while the extra- 
light and NSCs can overlap in morphological properties, 
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Fig. 19.— Mass vs. size for NSCs collected in Seth et al. (2008a). 
The large solid star shows the position of the NGC 404 NSC, show- 
ing that it has a size and mass typical of other early type galaxies. 
Data is from Carollo et al. (1997, 1998, 2002), Boker et al. (2002), 
Cote et al. (2006) and Seth et al. (2006). Points with underlying 
gray dots have dynamical or spectroscopic M/L estimates from 
(Watcher et al. 2005) and Rossa et al. (2006). 

the two can still be distinguished based on their kine- 
matics and stellar populations. We disagree with their 
argument. Kinematically, both NSCs and extra-light can 
be disky and rotating (Krajnovic et al. 2008; Seth et al. 
2008b; Trippe et al. 2008) and in early-type galaxies, the 
stellar populations of many NSCs are not significantly 
different from their underlying galaxy (Cote et al. 2006; 
Rossa et al. 2006). We therefore suggest that the distinc- 
tion between extra-light components and NSCs may be 
ambiguous in some cases such as NGC 404. 

The extra-light components in elliptical galaxies are 
thought to result from gas funneled to the center of the 
galaxy during gas-rich mergers (e.g. Mihos & Hernquist 
1994; Hopkins et al. 2009). Given the complicated kine- 
matics and stellar population of the NGC 404 NSC, it 
must be formed from multiple episodes of nuclear accre- 
tion. However, it formed about half its stars 1-3 Gyr 
ago. In del Rio et al. (2004), they suggest NGC 404 
underwent a merger with a gas- rich dwarf ~1 Gyr ago 
to explain the presence of the HI gas disk observed at 
large radii. The age estimate of the merger was based 
on the radius beyond which the gas had not yet settled 
onto the galaxy plane. It is therefore plausible that the 
~ 6 x 10 6 M© of stars formed at about this same age 
in the NSC resulted from accumulation of merger gas at 
the center of the galaxy. A merger origin is also consis- 
tent with the presence of counter-rotation in the central 
excess. This scenario is quite different from the forma- 
tion mechanism suggested by observations of late-type 
galaxies (Seth et al. 2006; Walcher et al. 2006; Seth et al. 
2008b), which appear to form primarily by episodic ac- 
cretion of material (gas or stars) from the galaxy disk. 
Resolved observations of a larger sample of NSCs in 
early-type galaxies is required to determine if NGC 404 
is an unusual case, or if NSCs in early-type galaxies typ- 
ically show some evidence gas accreted from mergers. 

7.2. The Possible NGC 404 MBH 

Our dynamical modeling in §6 provides mixed evidence 
for the presence of an MBH in NGC 404 with mass less 
than 10 6 M (making it an IMBH candidate). If verified 
through follow-up observations, this BH would be the 
lowest mass dynamically detected BH at the center of a 
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galaxy, with previous determinations of low mass BHs 
in M32 (2.5 ± 0.5 x 10 6 M Q ; Verolme et al. 2002) and 
Circinus (1.7 ± 0.3 x 10 6 M Q ; Greenhill et al. 2003) be- 
ing significantly more massive. A reverberation mapping 
mass of the BH in nearby Seyfert 1 galaxy NGC 4395 (Pe- 
terson et al. 2005), and indirect mass measurements in 
other low mass Seyfert 1 galaxies (e.g. Greene & Ho 2007) 
provide additional evidence for the presence of IMBHs 
at the centers of lower mass galaxies. Possible IMBHs at 
the center of star clusters (and possible former nuclei) Gl 
(Gebhardt et al. 2005), w Cen (Noyola et al. 2008; van 
dcr Marel & Anderson 2009), and M54 (Ibata et al. 2009) 
suggest that IMBHs with masses of 10 4-5 M may be 
widespread, but these results remain controversial. Mea- 
suring the occupation fraction and mass of BHs in galax- 
ies like NGC 404 is key to understanding the formation 
of MBHs in general (Volonteri et al. 2008). 

From the latest version of the Mbh — cr relations for 
different galaxy samples from Giiltekin et al. (2009), the 
predicted Mbh ranges from 0.8 — 3 x 10 5 M for a bulge 
cr of 35-40 kms -1 . This range of masses is consistent 
with all our dynamical BH mass estimates. A some- 
what higher mass (8 x 10 5 M ) is implied by the Haring 
& Rix (2004) relationship (see Fig. 18) and is consis- 
tent only with the gas-dynamical estimate, while the Fer- 
rarese et al. (2006) relationships have masses > 10 6 M Q 
and are thus inconsistent with our data. If a black hole 
with a few x 10 5 M resides in the center of NGC 404, 
better knowledge of the mass model from future HST 
observations should enable us to detect it. 

The relative mass of NSCs and BHs within a single 
galaxy can probe the relationship of these two objects. 
Seth et al. (2008a) find that the relative mass of BHs 
and NSCs is of order unity (with a scatter of at least 
one order of magnitude), while Graham & Spitler (2009) 
suggest a possible evolution in the relative mass with the 
least massive galaxies being NSC dominated, and the 
most massive galaxies being BH dominated. The BH in 
NGC 404 is <10% of the NSC mass and thus is not an 
obvious outlier from other galaxies where the mass of 
both objects is known. 

Finally, we note that the evidence of variability in the 
UV spectrum (Maoz et al. 2005) and possibly the NIR 
dust emission (§3.1), suggests that the NGC 404 black 
hole accretion is variable. In both cases, the observed 
flux decreased from observations taken in the 1990s to 
those taken in the last decade. This might help ex- 
plain the ambiguous nature of the AGN indicators in 
the NGC 404 nucleus; emission that originates near the 
BH (e.g. hard X-ray emission, broad Ha) may have re- 
cently disappeared, while the narrow-line region located 
at tens of parsecs, may continue to show signs of the ear- 
lier accretion. However, the observed hot dust emission 
suggests that the BH may still be accreting at a very low 
level. 

8. CONCLUSIONS AND FUTURE PROSPECTS 

This paper is the second resulting from our survey of 
nearby nuclear star clusters, and demonstrates the rich 
detail we can obtain for these objects. The NGC 404 nu- 
cleus is a complicated environment, with both a nuclear 
star cluster and a possible black hole. We have found 
that the surface brightness profile of the inner part of 
NGC 404 suggests the galaxy can be broken into three 



components: (1) a bulge that dominates the light be- 
yond 1", with M bu i ge - 9 x 10 8 M , r eff = 640 pc, and 
a Sersic index of ^2.5, (2) a NSC that dominates the 
light in the central arcsecond with r e // = 10 pc and a 
dynamical mass of 1.1 ± 0.2 x 10 7 M , and (3) a cen- 
tral excess at r < 0'.'2, composed of younger stars, dust 
emission, and perhaps AGN continuum. NIFS IFU spec- 
troscopy shows that the NSC has modest rotation along 
roughly the same axis as the HI gas at larger radii, while 
the central excess counter-rotates relative to the NSC. 
Furthermore, molecular gas traced by H2 emission shows 
rotation perpendicular to the stellar rotation. A stellar 
population analysis of optical spectra indicate that half 
of the stars in the NSC formed ~1 Gyr ago. Some an- 
cient and very young (<10 Myr) stars are also present. 
This star formation history is dramatically different from 
the rest of NGC 404, which is dominated by stellar pop- 
ulations >5 Gyr in age. We suggest a possible scenario 
where the burst of star formation in the NSC ~1 Gyr 
ago resulted from the accretion of gas into the galaxy 
center during a merger. This formation scenario is quite 
different from the episodic disk accretion suggested by 
observations NSCs in late-type galaxies. 

Our dynamical modeling of the stellar and gas kine- 
matics provide mixed evidence for the presence of a 
black hole in NGC 404. Assuming a constant M/L 
within the nucleus, the stellar dynamical model suggests 
an upper limit of 1 x 10 5 M , as well as measuring a 
M/L/ = 0.70 ±0.04 for the NSC. The gas kinematics are 
best fit by models including the presence of a black hole 
with M B h = 4. 5^0 x 10 5 M . Both dynamical BH mass 
estimates rely on a model for the stellar mass that we 
construct from HST F8141F-band imaging. Uncertain- 
ties in the light profile (due to variability) and M/L (due 
to stellar population changes) within the central excess 
are of the same order as the difference between the two 
black hole mass estimates. We have proposed to mea- 
sure the mass model by using HST/STIS spectroscopy 
to resolve the stellar populations within the nucleus and 
additional multi-band imaging to extend this model to 
two dimensions. If successful we will combine this mass 
model with the kinematic observations presented here as 
well as larger-scale kinematics obtained from the MMT 
to more robustly determine the BH mass. 

In addition to the direct evidence of the BH from the 
dynamical models, we find two other properties of the 
nucleus which suggest the presence of an AGN. First, 
we find unresolved hot dust emission at the center of the 
NSC with a luminosity of <~ 1.4 x 10 38 ergs s _1 . Compar- 
ison of the NIFS light profile to previous NICMOS obser- 
vations suggests that this dust emission may be variable 
as is seen in other AGN. Second, the H2 line ratios within 
the central arcsecond indicate thermal excitation in dense 
gas similar to what is seen in other AGN. Our proposed 
HST observations request multi-epoch UV through NIR 
imaging to search for definitive evidence of BH accretion 
in NGC 404. 

Our nearby NSC survey includes 13 galaxies within 
5 Mpc with Mb between -15.9 and -18.8 for which we will 
have comparable data to that presented here obtained 
using MMT, VLT, and Gemini. These galaxies span a 
wide-range of Hubble types in which we can examine the 
process of NSC formation; there are four early type E/S0 
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galaxies and 9 spirals of type Sc and later. NGC 404 
represents one of the stronger cases for finding a black 
hole, given its LINER emission and proximity. However, 
we expect to be able detect or place upper limits of < 
10 5 M© on black holes in each of the sample galaxies. 
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